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Abstract 

In this paper, we consider a stabilization problem of an uncertain system in a networked control setting. Due to the network, the 
measurements are quantized to finite-bit signals and may be randomly lost in the communication. We study uncertain autoregressive 
systems whose state and input parameters vary within given intervals. We derive conditions for making the plant output to be mean square 
stable, characterizing limitations on data rate, packet loss probabilities, and magnitudes of uncertainty. It is shown that a specific class of 
nonuniform quantizers can achieve stability with a lower data rate compared with the common uniform one. 


1 Introduction 

This paper studies stabilization of a linear system in which 
the plant outputs are transmitted to the controller through a 
bandwidth limited lossy channel and the exact plant model 
is unavailable. For control over finite data rate channels, it is 
well known [1,2] that there exists a tight bound on the data 
rate for stabilization of linear systems, which is expressed 
simply by the product of the unstable poles of the plant. 
Such data rate limitations have been developed under a va¬ 
riety of networked control problems. For general nonlinear 
systems, it has been pointed out that the limitation is related 
to topological entropy [3-5]. For an overview on the topic, 
we refer to [6]; for more recent works, see, e.g., [7,8]. On 
the other hand, control over packet dropping channels has 
also been studied actively (see, e.g., [9,10]). Interestingly, 
by modeling the behavior of the losses as i.i.d. random pro¬ 
cesses, the maximum packet loss probability for achieving 
stabilization can also be characterized solely by the product 
of the unstable poles of the plant. Recent works have ex¬ 
tended such results to the case of Markovian packet losses. 
In particular, [11] and [12] have derived the minimum data 
rates for the static and time-varying rate cases, respectively. 
We note that the works mentioned above assume perfect 
knowledge of the plant models. 
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Despite the active research in the area of networked con¬ 
trol, uncertainties in plant models have received limited at¬ 
tention and are thus the focus of this work. In general, it is 
difficult to deal with the combination of uncertainties in the 
systems and incompleteness in the communication. In par¬ 
ticular, in data rate limited control problems, the state evolu¬ 
tion must be estimated through quantized information, but in 
the uncertain case, this task becomes complicated and often 
conservative. In [13,14], linear time-invariant systems with 
norm bounded uncertainties are considered, and controllers 
to robustly stabilize the systems are proposed. In [15], scalar 
nonlinear systems with stochastic uncertainties and distur¬ 
bances are studied, and a data rate bound sufficient for the 
moment stability is derived. Moreover, related stabilization 
problems are studied from the viewpoints of adaptive con¬ 
trol [16] and switching control [17] as well. In these results, 
however, only sufficient conditions on data rates have been 
obtained, and they are not concerned with characterizing 
the minimum. On the other hand, observation problems of 
nonlinear time-varying uncertain systems are studied in [4]. 
Both required and sufficient data rates for observability are 
characterized by using the notion of topological entropy. 

More specifically, we consider the stabilization of a para¬ 
metrically uncertain plant over a Markovian lossy channel. 
The plant is represented as an autoregressive system whose 
parameters vary within given intervals. We develop bounds 
on the data rate, the packet loss probability, and plant un¬ 
certainty for stabilizability. The results become tight for the 
scalar plants case. In the course of our analysis, we demon¬ 
strate that the data rate can be minimized by employing a 
class of nonuniform quantizers, which is constructed in an 
explicit form. These quantizers have an interesting property 





that the cells are coarser around the origin and finer further 
away from the origin; such a structure is in contrast to the 
well-known logarithmic quantizer [18]. 

For the case of uncertain state parameters, the authors have 
studied the minimum data rate under model uncertainties and 
packet losses for the i.i.d. [19] and the Markovian [20] lossy 
channels. It is also noted that, in [21], for a similar class of 
uncertain plants, stabilization techniques have been devel¬ 
oped based on the logarithmic quantizers [18]. This paper 
aims at further studying the more realistic situation where 
uncertainty is also present in the actuator of the plant. That 
is, the parameter of the control input may also be uncertain. 

The main difficulty in the current setup can be described 
as follows. Evaluation of the estimation error and its evo¬ 
lution is the key to derive the minimum data rate. Due to 
plant instability, the estimation error grows over time, but it 
can be reduced based on state observations. In our previous 
work [20], we have assumed uncertainty only in the state 
coefficients. In the presence of uncertainty in the actuator, 
the results there are not applicable. In particular, when the 
control input is large, the estimation error will grow further, 
making the analysis more involved. We show that uncer¬ 
tainty in the actuator side introduces additional nonunifor¬ 
mity in the quantizer structure when compared to our pre¬ 
vious results. 

This paper is organized as follows. In the next section, we 
formulate the stabilization problem for the networked con¬ 
trol system. In Section 3, we consider the fundamental case 
for the scalar plant systems. The general order plants case 
is considered and a sufficient condition for the stability is 
shown in Section 4. In [20], a necessary condition for the 
general order plants case is also provided. However, in this 
paper, we do not include the corresponding result since in 
general it contains some conservativeness and its signifi¬ 
cance may be limited. Finally, concluding remarks are given 
in Section 5. The material of this paper was presented in [22] 
in a preliminary form, but this version contains updated re¬ 
sults with their full proofs. 

Notations: Z+ is the set of nonnegative integers and N stands 
for the set of natural numbers. Iog 2(0 is simply written as 
log(-). For a given interval J/ on R, denote its infimum, 

supremum, and midpoint by J/, J/, and c(J/) + J/)/2, 

respectively; its width is given by ;u(J/) 

2 Problem setup 

We consider the networked system depicted in Fig. 1, where 
the plant is connected with the controller by the communi¬ 
cation channel. At time k e Z+, the encoder observes the 
plant output yt e R and quantizes it to a discrete value. 
The quantized signal s^ e is transmitted to the decoder 
through the channel. Here, the set Zff represents all possible 
outputs of the encoder and contains N symbols. Thus, the 



Fig. 1. Networked control system 

data rate R of the channel is given as R log A; thus, the 
term refers to the bit rate of the communication, which is 
consistent with the literature [11,12]. The decoder receives 
the symbol and decodes it into the interval ifk c R, which 
is an estimate of y^. The transmitted signal Sk may be lost in 
the channel. The result of the communication is notified to 
the encoder by the acknowledgment signals before the next 
communication starts. Finally, using the past and current es¬ 
timates, the controller provides the control input Uk e R. 

The plant is the following autoregressive system with un¬ 
certain parameters: 

yk+i^aikyk+a2,kyk-i + - ■ ■+an,kyk-n+i + btut. ( 1 ) 

The parameters are bounded and may be time varying as 

a^k G := [fl* - ^i, a* + ^i], i - 1,2,... ,n, 
bkeS:^[b* -6,b* + 6], (2) 

where e, > 0 and 5 > 0. The initial values yk, k - -n + 
1,..., -1,0, are in the known intervals as yk e Yk, where 
0 < /i(Ti) < oo. To ensure controllability at all times, we 
introduce the following assumption: For every time A: 6 Z+, 
the input parameter bk is nonzero. That is, 

\b*\ -6>0. (3) 


In the communication channel, transmitted signals may be 
randomly lost. Denote the channel state at time k by the 
Markovian random variable jk e (0,1]. This state represents 
whether the packet is received {jk - 1) or lost {jk - 0). 
The loss probability at time k depends on the previous state 
jk-i and is denoted by the failure probability p and the 
recovery probability q as follows: Prob(y,t = 01 jk-i = 0) = 
1 - q, Prob(7,t = 1 | jk-i = 0) = ProbCy^t = 0 | jk-i = 1) = 
p, Prob( 7 i: = 1 \ yk-i = 1) = 1 - /?. To make the process 
{yk}k ergodic, assume p,q € (0,1). Moreover, without loss 
of generality, assume that at the initial time k - 0 the packet 
is successfully transmitted, i.e., yo = 1- 

The communicated signal Sk e Eat := (1,2,..., A] is the 
quantized value of the plant output yk and is generated by 
the encoder as Sk - (f>N(yklo'k)- Here, the quantizer 0v(') is 
a time-invariant map from [-1 /2,1 /2] to Eat with a scaling 
parameter crk > 0. The quantizer divides its input range 
[-1/2,1/2] into A cells and its output is the index of the cell 
into which the input falls. We assume the boundaries of the 
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quantization cells to be symmetric about the origin. When 
N is even, we denote the boundary points of nonnegative 
quantization cells as hi, I - 0,1,, [N/T], where 

ho = 0, /ifw/21 = 2’ ^ 

If N is odd, the nonnegative boundaries can be written by hi, 
I - 1 , 2 ,..., INjlX where 0 < hi < h 2 < ■ ■ ■ < /jfA^/ 2 i = 1/2. 
For notational simplicity, we add ho - 0 and use the same 
notation for this case also. 

Based on the channel output ytSk, the decoder determines 
the interval c R, which is the estimation set of y^. When 
the packet arrives successfully, i.e., yk - 1, is the quan¬ 
tization cell which falls in. Otherwise, is taken as the 
entire input range [-cr^/2, cr^/2] of the quantizer. 

The initial value of the scaling parameter cr^. and its update 
law are shared between the encoder and the decoder. The 
parameter cr^. is updated as follows. At time k, the encoder 
and the decoder predict the next plant output y^+i based on 
the past estimates ^k-i+i, i - 1,2,..., n. As time progresses 
from klok+\, the past outputs yt-M e d/jt-i+i are multiplied 
by ai^k G - Let b® the set of predicted values for 

yt+i in the form of cii^tyk-i+i and is given by 

^k+i {‘^13'i: + • • • + any'k^n+\ • e e ^n, 

^ ^k, ■ ■ ■ ,yk-.n+\ ^ ^k-n+\}- (5) 

Moreover, since the applied input is bkUk, the set {y” + but : 
W G }f^^i,b e S) is large enough to include y^+i. Note that 
the above prediction set in (5) is computable on both sides 
of the channel by the acknowledgment signal regarding y^-i 
from the decoder to the encoder. Finally, this set {y“ -H buk : 
y^ € Z? e Sj must be covered by the quantizer to avoid 
saturation. Hence, the scaling parameter cTk is a function of 
and must be large enough that 


the scaling parameter crk satisfying (6) and a controller (7) 
such that the worst case output yk over all deterministic 
perturbations is mean square stable (MSS): E[sup^gj/^ |y|^] ^ 
0 as —> oo. Here, ^k is the decoder output at time k and 
the expectation is taken with respect to the packet losses yo, 

■ ■■,yk. 

Remark 2. We note that the results and the proofs in this 
paper have been polished compared with our previous paper 
[20] thanks to the help of the anonymous reviewers. In [20], 
the definition of stabilizability and the proof of the necessary 
condition for the scalar plants case should be updated. This 
however does not change the bounds on the data rate and 
the loss probabilities shown in the theorem. 

Remark 3. To simplify the analysis in dealing with un¬ 
certainties, we have introduced some structures in the en¬ 
coder and the controller. While they are closely related to 
those employed in, e.g., [13,15,23] for obtaining sufficient 
conditions, there is some conservatism. In Definition 1, the 
supremum is taken over ifk, which contains all possible yk 
over aij e f = 1,2,..., n, bj e S, and 

yj e Yj, j - -n + 1, -n H- 2,..., 0. If we do not limit the 
controller class to (7), there may exist one which can com¬ 
pute an estimation set tighter than J/^, though it is difficult 
to describe the tightest estimation set analytically. Further¬ 
more, for Definition 1, it is important that the quantizer does 
not saturate. This is guaranteed by (6), and there is always 
a quantization cell containing yk- 

3 Scalar plants case 

We first analyze the simple setup with the scalar plant (a 
first-order autoregressive process): 

yk+i = akyk + bkUk, 

akeJl:=[a-e,a+e], bk e S ^ [b*-S,b'‘+6], (8) 


o-k+i > 2 sup |y H- buk\. (6) 

beS 


The controller provides the control input ut based on the 
past and current estimates ^k-n+i ,..., as 

n 

M/t = ^ .fi.k , (7) 

1=1 

where /i,i(-) are maps from an interval in R to a real number, 
which determine the input based on the estimates. 

This paper investigates stabilization of the uncertain net¬ 
worked system in Fig. 1 by designing the encoder, the de¬ 
coder, and the controller under the constraints (4)-(7). 

Definition 1. The feedback system depicted in Fig. 1 is 
stabilizable if there exists a pair of an encoder <pN with 


where e > 0 and 6 > 0. We assume that its dynamics is 
always unstable in the sense that the parameter Ok has mag¬ 
nitude greater than 1 at all times, i.e., 

|fl*|-e>l. (9) 


To express the main result of this section, let 

\b*\-6 


\a*\ + e 


n := 


, A:=e-H(5—, 

\b*\ + 6 \b*\’ 


V -t 1 + 


/7{(|a1 + e)2-l) 


l-(l-^)(|a1 + e)2' 


( 10 ) 


Here, r^, r^, and A reflect the magnitudes of the uncertainties 
and V represents the effect of packet losses in the required 
data rate as we will see in (11). We show later, in the proof 
of the next theorem, that the radicand of v is positive when 
the feedback system is stabilizable. If ak is a constant, i.e.. 
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Ok - a for all time k, it has been shown in [11] that mean- 
square stability implies that ^ > 1 - 1/a^. Taking account 
of the case that a -\a*\ + e, we have > 0. The following 
theorem shows a condition on the data rate R - log N, the 
loss probabilities p, q, and the magnitude of uncertainty A 
for stabilizability. 


Theorem 4. Consider the feedback system in Fig. 1 with 
the scalar plant in (8). If the system is stabilizable, then the 
following inequalities hold: 


R ^ Rnec ■ — 


^ ^ ^nec 1 ■ 


log 

(log|a*| -Hlogv 

1 


(|fl1 + e)2 




if e>0 or 6>0, 
if e - 6 - 0, 

p^{l-(\a*\ + er^} 

1 - A2 


0 < A < 1. 


( 11 ) 

( 12 ) 

(13) 


Furthermore, if these inequalities are satisfied with an even 

N, it is possible to construct a stabilizing controller. 

We note that when (12) and (13) hold and e > 0 or d > 

O, we have that 0 < 1 - Av < 1 and hence /?nec is well 
defined; see the last part of the proof of Theorem 4 for 
details. This theorem provides limitations for stabilization 
on the data rate, the packet loss probabilities, and the plant 
uncertainty. The required data rate /?nec and the recovery 
probability q^c are monotonically increasing with respect to 
the uncertainty bounds e and 6. This means that more plant 
uncertainty requires better communication with higher data 
rate and recovery probability. 


We see that the sum of the uncertainties A = e H- 6\a*\l\b*\ 
appears in the limitations. It is interesting that there is no 
explicit limitation on e or 6, but the sum A of these uncer¬ 
tainties must be smaller than 1. This indicates some trade¬ 
off in the tolerable uncertainties for at and bt. In particular, 
the product d|a*| implies that for more unstable plants, the 
bound 5 on the input parameter bt has more effect on the 
stability conditions in the theorem. We remark that when 
the input parameter is known and is constant as bt = b*, 
i.e., d = 0, then the limitations /?nec and coincide with 
those shown in [20], where uncertainty is present only in 
at. Moreover, if at is also known, i.e., e = <5 = 0, then the 
limitations are equal to those in [11], where the exact plant 
model is assumed to be available. 


We provide an example to illustrate the limitations in The¬ 
orem 4. Consider a plant with a* - 2.0 and b* - \.Q and a 
channel with the loss probabilities p - 0.05 and q - 0.90. 
Fig. 2 shows the bound /?nec on the data rate versus the un¬ 
certainties e in at and 6 in bt. When the sum of the uncer¬ 
tainties A is large as (12) is not satisfied, 1 - Av in (11) is 
nonpositive and hence the required data rate for stabilizabil¬ 
ity becomes infinite. 



0.1 0.2 0.3 0.4 0.5 0.6 

Magnitude of uncertainty e in at 


Fig. 2. The data rate limitation R^ec versus the magnitudes of the 
uncertainties e and S (a" = 2.0, b" = l.O, p = 0.05, and q - 0.90) 


nite data rate and lossless channels. We remark that those 
conditions contain conservatism even for the scalar plants 
case. Consider the scalar plant (8) when the input coefficient 
is known (d = 0) and the channel is lossless. For such sys¬ 
tems, the sufficient bound on the data rate R^ut in [13] and 
the one from [15], respectively, become 


^suf log 


\a*\ - e(|a*| -H e) 

1 - e(2|fl*| + 2e+l) 


and ^Lf-log 


It can be verified that our result is tighter than these bounds 
as /?nec < ^suf and /?nec < For general order plants, how¬ 
ever, it is difficult to compare these results since the types of 
uncertainties are different: In [13], unstructured uncertain¬ 
ties are considered, and it is hard to describe the data rate 
limitation in an explicit form, while [15] deals with nonlin¬ 
ear plants but only scalar ones. 


We now present the proof of Theorem 4. The key idea lies in 
evaluating the expansion rate of the state estimation sets due 
to plant instability. The proof consists of two steps, which 
are presented in Sections 3.1 and 3.2, respectively. 


Remark 6. Compared with our previous work [20], the main 
difficulty is that we have to take account of the expansion in 
the state estimation sets by control inputs. If we know the 
exact control input applied to the plant, then the width of the 
estimation set is not affected by the input since we can track 
the variation of the state precisely. However, in the current 
setup, the estimation set may expand by the control input 
due to the uncertainty in bt. Hence, the scaling parameter 
crt+\ must be selected to cover this expansion in addition to 
that by plant instability. 


3.1 The quantizer minimizing the expansion rate 


Remark 5. The works of [13] and [15] have shown suffi¬ 
cient conditions for stabilization of uncertain plants via fi¬ 
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In this subsection, we introduce the expansion rate for a 
given quantizer. Then we show the optimal quantizer which 
minimizes the rate in the worst case. 



For a given quantizer whose boundary points are let 


i e such that hi - hi, h^ - h^, and w/ are the same 

for all / e {/, / + 1,..., m - 1}, i.e., for a constant w 


Wl 


2(|a*| + e)/j;+i if N is odd and I — 0, 

< (\a*\ + e){l + jir^hM-i\a*\-€){l-jlr^h, 
else, 


(14) 


Wl - w, VI e [i,i + I,... ,m - 1], (19) 

We shall show that for all quantizers it holds that 


for 1 = 0,1,, riV/2] - 1. We see later in Section 3.2 that 
this Wl characterizes the expansion rate of the volume of the 
estimation set for one sampling period due to the uncertain 
parameters at & ^ and bt e S. The rate varies depending 
on I, which represents the cell which observed output falls 
into and we have to consider the worst case to guarantee 
stability against uncertainties. 


max wi(h) > w, (20) 

where wi(h) denotes the expansion rate w/ in (14) with the 
quantization boundaries {hi]’ll.. This is done by contradic¬ 
tion. Suppose that max;e(;;+i m-uwi(h) < w. Then, from 
(19), it follows for all I e {/, /'+' 1,..., m - 1) that 


The tight lower bound on the worst-case expansion rate 
can be derived as shown in Lemma 7 below. Let us de¬ 
fine the boundary points dividing [-hm,hm], where 

m e {1,2,..., [Af/2]), as 


K := 


(u i-t(rarb)' 
i \-t(ranf 

I 

\ m—t' 


if e > 0 or (5 > 0, 
if e = d = 0. 


(15) 


Furthermore, define w^, which is used in the following 
lemma to represent the worst-case expansion rate as 


w 


m * 


t := 


t 


|/i„(|fl*| + e)(l + ^) 

"t m—t' 




if N is odd, 

1 if is even. 


i if is odd, 
0 if Af is even. 


if e > 0 or d > 0, 
if e = d = 0, 


(16) 


wi(h) < max wi'(h) < w - wi(h). (21) 

We shall compare hi with hi for each / using (21). For I - i, 
we have hi = /?,. Substituting these into (14) yields 


2(\a*\ + e)hi+i if N is odd and i - 0, 

Wi{h) = < (|a*| -H e) (l -H hi+i - (|a*| - e) (l - j^) hi 

else, 

1 2(\a*\ + e)hi+\ if N is odd and i - 0, 

(|fl*| + e) (l + ^) hM - (|a*| - e) (l - ^) /),■ 

else. 


Since Wi{h) < Wi(h) from (21), it follows that 

hM<hi+i. (22) 

Moreover, by (14), we have for / = / -H 1, / H- 2 ..., m - 1 that 


Lemma 7. Given a quantizer dividing [—1/2,1/2] 

into N cells, consider a subset of the quantization region 
[—hm,hm], where m e (1,2,..., [A/^/2]}. The worst-case ex¬ 
pansion rate of the cells in [—h^,h^] is bounded as 


wi{h) = i\a*\-ve)^l-v-^^hi+i - (|a*|-e)^l-|^j/i;. 


wi(h) = (|a*|-i-e) 


max 


Wl > wl,. 


(17) By substituting these into (21), we obtain 


The equality in (17) holds if hi — h^ for I — 0,1,..., m. 
Furthermore, consider the subset of the quantization region 
[hm',hm], where I < m' < m < \NI2\ Then, it follows that 


max Wl 

/G(m', m' + l, 1 } 


{/*«-} (|fl*| + e) (1 + 


ife^d^Q. 


-rgrt 

if e > 0 or 6 > 0, 


(18) 


Proof. To prove the inequalities (17) and (18), we first 
assume that there exists a set of boundary points 1^;}™,-, 


hi+i < rarbhi -H 


h+\ = Vbhi + 


max/-e|i.i+i,...,m-i| wr(h) 

(\a*\ -H e)(l -H j^) 
w 

(|fl*|-He)(l -I- j^) 


By introducing the relation (22) to the above, we recursively 
obtain hi < hi for every I e {i-vl,i-\-2,..., m]. This however 
is in contradiction with h^ — hm- Thus (20) holds. 


To establish (17) and (18), we must show the existence of 
{ft;}™. satisfying the assumptions made at the beginning of 
this proof. First, for (17), consider the quantizer {ftJl^Q in 
(15) dividing the subset [-ftm, ftm] of the quantization range. 
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Fig. 3. Boundaries of the optimal quantizer (N = 8, a* = 3.0, 
e = 0.5, andZ?* = 1.0) 



(a) Uniform quantizer case (b) Optimal quantizer case 


Fig. 4. Expansion of quantization cells: With the uniform quantizer, 
cells further away from the origin expand more at each time after 
applying control, while with the optimal one, all cells result in the 
same width. 


By a routine calculation and the relation /to = 0 in (4), one 
can confirm that {/r*)™o satisfies - ho - 0, /t*, = h^, and 
(19) with i = 0 and tv = tv*,. Hence, (17) follows with the 
equality condition. 


Next, to show the case of (18) with [/tms /tm], suppose that i — 
m' and consider the following sequence for I - 0,1,... ,m- 
m': 


^m'+l — * 


{rarb)‘hm' + [hm - (farbT ’"'hm'] 


hfn' "h (.hfyi hffi' ') 


I 


if e > 0 or d > 0, 
if e = (5 = 0. 


We then have that h^’ - h,„>, - h^, and (19) where tv 

is equal to the right-hand side of (18). This concludes the 
proof. □ 


Let us denote by the quantizer consisting of the bound¬ 
aries in (15). This quantizer is optimal in the sense 

that it minimizes maXie|o,i,...,fA?/ 2 i-i) w;, and the minimum is 
given as from (17). Quantization of is nonuniform 
and becomes more so when the plant has more uncertainty. 
To see this, consider the plant with a* - 3.0, e = 0.5, and 
b* - 1.0, and take = 8. The boundaries of for 6 = 0.0 
are shown in Fig. 3 (a) and for 6 - 0.3 in Fig. 3 (b). 

The nonuniformity of (p*^ is an outcome of minimizing the 
effect of the plant uncertainty on the state estimation. This 
characteristic can be explained as follows. Due to quanti¬ 
zation, only the interval J/j, containing the true output is 
known to the controller. After one time step, because of the 
plant instability, the interval in which the output should be 
included will expand. When the plant model is known, the 
expansion ratio is constant and is equal to |a*| for any quan¬ 
tization cell. However, with plant uncertainties, the ratio de¬ 
pends on the location of the cell. In particular, cells further 
away from the origin expand more. Moreover, to bring such 
cells around the origin, larger control inputs are required 
compared with cells closer to the origin. Because of the un¬ 
certainty in the input parameter bk, larger control input will 
result in further expansion of the interval. This fact is illus¬ 
trated in Fig. 4(a) when uniform quantization is employed. 
In contrast, when the proposed quantizer 0^ is used, the cells 
after one step have the same width (Fig. 4(b)). We note that 


when there is no uncertainty in the plant, i.e., e = d = 0, 
then becomes the uniform quantizer. 

The quantizer structure above is similar to that introduced 
in our previous work [20], which studied plant uncertainties 
only in aiy. We also note that nonuniform quantizers have 
appeared in the networked control literature, e.g., in stabi¬ 
lization problems [18,24] and identification problems [25]. 

3.2 Derivation of the limitations in Theorem 4 

We are now ready to prove Theorem 4. The central problem 
in the proof is to evaluate the expansion rate of the estima¬ 
tion set during one sampling period. First, a tight bound 
on the expansion rate is shown, and then with the bound 
and the optimal quantizer the limitations (11)-(13) for 
stabilizability are derived. 

Proof of Theorem 4. (Necessity) Suppose that the feedback 
system is stable in the sense of Definition 1 with an encoder 
and a controller. We first show that E[sup^gj^^ |yp] ^ 0 as 
A: —> oo implies that E[cr^] ^ 0. The estimation set 3fk avail¬ 
able at the controller corresponds to a quantization cell. By 
(6), this set ffk is guaranteed to contain jk. Let d > 0 de¬ 
note the smallest width of the quantization cells. Then, from 
the definition of the quantizer, we have that ^i{3fk) ^ dcrk, 
and hence sup^.^^/^ |j| > dak 12. Therefore, lim^^oo E[cr^] = 0 
follows from the mean square stability of the system. 

The rest of the proof consists of three steps. The first step is 
to prove the following inequality, which provides a bound 
on the expansion rate of cr^: 


CTjt+l > TJkak. 


Here, rjk is the random variable defined as 

j\a*\ + e if 7 ^ = 0, 

m := 1 -f 1 

[wi, ifn = i. 


(23) 


(24) 


where w/ is defined in (14) and /,i e [0,1, • ■ ■, \N/2] - 1} is 
the integer such that infj,£j/j \ylak\ - hi,., which is the index 
of the quantization cell which yk falls into. 
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To establish (23), we recall ( 6 ) and claim that for any control 
input Uk, we have that 

cTk+\ >2 sup |y“ + bukl 

bes 

>2 sup |y“+ 

beS 

= (25) 

where the input m*(-) is defined as 


infimum as = {a* + e)^k and = (a* - e)^k, 

respectively. Substitution of these into the left-hand side of 
(28) gives 

+ jp| lc(J/^+i)l 

= {(a*-He)(l-H^)Vi -(a*-e)(l-p)/z4|cr^. 
Hence, by (14) the relation in (28) holds. 


U0) 


c(J/) 

ft* 


(26) 


for an interval J/ on R. The input m*(J^^^_i) brings the mid¬ 
point of the prediction set into the origin when the 
parameter bt is equal to b*. 


For the derivation of (25), we first consider the case that 
c(J/^^i) > 0 and b*-S > 0. In this case, it is obvious that the 
input Uk minimizing supj,-gj/- \y^ + buk\ is nonpositive. 
For such inputs ut < 0, we have that 


sup ly + buk\ 

= max {^,^1 + {b* - 6)uk, -^kyx “ (^* + • (27) 


(ii) J/it < 0 < J/it: In this case, we have 

= (a* + e)^,, = (a* + e)^^. (29) 

Moreover, because of the symmetry in the quantization 
cells, it holds that - 0 and hence m*(2^^^i) = 0. 

Thus, cr*^j = To compute this width /i(J/^^j), 

consider the following two cases: (ii-1) If jk - 0, 
then = [-crj./2, cr^/2]. Hence, by (29), we have 
/i(J/^^i) = {a* + e)crk- (ii-2) Otherwise, N must be odd and 
4 = 0 from the symmetry of the quantizer and the condition 

(ii) . Thus, ;u(d/^^j) = 2(a* + e)h\crk- Hence, (28) holds for 
this case also. 

(iii) < 0: This case can be reduced to (i). 


Take an arbitrary nonpositive input and denote it as Uk - 
M*(J/^^[) -I- A„j. With this expression and (26), the right-hand 
side of (27) is equal to 


2 




2 


^c(J/,-,,) + (f.*-d)A„„ 
^c(J/,-,,)-(f.* + 5)A„,|, 


which takes its minimum when A„j = 0. This proves (25). 
The case of c(J/^^j) < 0 or Z?* -15 < 0 can be reduced to the 

above by appropriately flipping signs of and b*. 

The lower bound cr^^j on (Jk+\ in (25) is in fact equal to the 
right-hand side of (23), i.e., 


From (25) and (28), we have established (23). 

As the second step, we consider the mean squares of both 
sides of (23) and derive an inequality regarding the expan¬ 
sion rate of crk- Note that due to packet losses, the expan¬ 
sion rate rjk is a random variable that depends on the previ¬ 
ous channel state. This time dependency causes difficulties 
in the analysis. To avoid this, we consider intervals between 
the times at which packets successfully arrive at the con¬ 
troller side; it is known that the intervals become an i.i.d. 
process [27]. We formally state this fact in the following. 

Let tj, j e Z+, be the times at which packets arrive, i.e., = 

1. From the assumption on the initial communication, stated 
in Section 2, we have 0 = to < 4 < • • ■ • Then, denote the 
interval between the times tj_i andt^byzy 4 - 4 - 1 , j > 1 . 
The process {Tj]j is i.i.d. and, for each j, it holds that 


^k+i = ^ko-k- (28) 

Here, tjk is present because the expansion rate of crk is af¬ 
fected by the packet losses yk- To derive (28), we consider 
three cases (i)-(iii) depending on the location of the estima¬ 
tion set ^k as follows. For simplicity, assume a*,Z>* > 0. 

(i) ifk > 0: This case occurs only when the packet arrives, 
i.e., jk - 1. In addition, from (24), rjk - wi^ where N is even 
or Ik + 0. From the basic results for products of intervals [26] 
and (9), for the interval we obtain its supremum and 


Prob(Ty = i) - 


l-p 

pq( l - qT^ 


if / = 1 , 
if / > 1 . 


(30) 


From (23) and the fact that 


n 


1 if A: = 4 , 

0 if A: G [ty + 1 , ty'+l), 


a lower bound on the expansion of crk from time 4 to 4+1 
is given by > (\a*\ + Here, the value of 
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w;,, varies depending on l,j, which is the index of the cell 
that ytj fell in. Since we have to take account of all possible 
deterministic perturbations, we consider the maximum over 
the suffix of wi,^: 


o-t .^1 > (|a*| + ' max WiO-,., (31) 


where contains the set of all possible I,., which is a subset 
of the index set {0,1,, riV/2] - 1). More precisely, at time 
k, Lk is the indices of the quantization cells which intersect 
with {y“ + buk-\ : e b e S). 


As to the right-hand side of (31), we shall prove that 


case. When e > 0 or <5 > 0, a routine calculation shows that 


1 


sup |y I - (r^ri)"* inf \y \ \ - 

\y-eif- \i-(rartr 


> 


2 ^ \b*\ 


Y 


1 


■ tirari,)" 


2{1 - tirarb)"'] 


(34) 


Here, for the first inequality, we fiave used sup^,_g^_ \y \ - 

^i0,)|2 + \c0;)\, inf^,_,^^ |y-| = -^i0;)/2 + |c(ifY)|, and 

0 < tr < 1, and the second one follows by 2 J/^. From 
(33) and (34), and by m < rA^/21, we arrive at (32). 


maxwio-k > Wp^/2-|0-^ (32) 


where w* and cr^ are given by (16) and (25), respectively. 
Let be the union of the quantization cells which have 
intersection with (y“ + buk-i : y“ e ^ Consider 

the following two cases. 


From (25), (31), and (32), we have 

> E[(|fl*| + (wf^/ 21 )" E[</]. (35) 

Notice that Tj+i is independent of cTty Since E[cr^] ^ 0 and 
by (25), the coefficient of E[cr* ^] in (35) is smaller than 1. 
After some computation using (30) and that {Tj]j is i.i.d., 
we obtain 


(i) vnf^k ^ 0 ^ supdf^: Let m e {1,2,..., [A/2]} be 
the cardinality of Lk. We then have that max/g^t wi — 
max/e|o,i,...,m-i| wi and sup^^^^- |y“| = h,nO-k- Thus, from (17) 
in Lemma 7, it holds that 


maxwia-k >-jf- sup |y | > 2 wpjv/ 2 i Id I’ 

y-'^, 


E[(|a*| H- e)2T>+i-0] = 1 h- 


p[{\a*\ + ef - 1 } 

1 - (1 - q)(\a*\ + ef 


and 1(1 - q){\a*\ + eY\ < 1. The inequality boils down to 


0 < (1 - ^)(|a*| H-e )2 < 1 , 

which implies > 0. Hence, we arrive at 


(36) 


where we have used > 2 wp^^ 2 .| to obtain the second 

inequality. Eurthermore, since 2 {y“ -H buk-i : e 

J/^, b e S] and by (25), we have 2 sup^_g^_ |y“| > cr*. Thus, 
we obtain (32) for this case. 

(ii) inf j/Y > 0 or supd/^ < 0: In this case, we can de¬ 
fine integers m' and m such that 1 < m' < m < rA/2], 
ly^l = h^’O-k, and sup^^^^- |y"| = h^crk- By using 
(18) in Lemma 7, it follows that 


ywf ^/21 < 1. (37) 

The last step is to derive (11)-(13) from (36) and (37). Eirst, 
suppose that e or d is positive and that N is even. Then, with 

the definition ( 16 ), it follows that Wp^^ 2 ] - 

Substituting this into (37), we obtain (jarh)^^'^ < 1 - vA and 
by taking logarithm of both sides, 

\og{rarh) 


max wiCFk 

leLk 

'{sup^_,^- ly-|-(/-„r,)'«-'inf^_,^- |y-|} 
x(|«-| + e)(l + j^)^-i-^ ife>0orh>0, 

ife = d = 0. 

\j^k' m-m' 


(33) 


Note that 


1 - vA > 0 


(38) 


since (farb)^^^ > 0. If A is odd, then similarly, from (37), 
we have that 


A > A<°> := 2 


log(l - Av) - logt 


- 1 , 


If e = d = 0, noticing that piS/k ) ^ we have (32) for this 


log(r„ri) 



where t (1 + 6l\b*\) j (\ - e/|a*|). Comparing the lower 
bounds and N^°\ by the assumption (3) on b* and (9), we 
can show N^°'' > N^'^\ Thus, the smaller bound appears 
as the data rate limitation in (11). To establish (12) and (13), 
we see from (36) and (38) that 

(1 - a 2) {l - (1 - q)(\a*\ + ef) > {(|fl*| + ef - l}. 


Here, w, is defined for the given quantizer as follows: 


max{wf\ wP) 

if N is odd and J?l, ^ 0, 

f , _(0)t 

max e,+d]^, w] 

if N is odd and J?l, 3 0, 

—(1) 
w: ' 

if N is even and $ 0, 

^ aKI 

-H djpi 

if N is even and 3 0, 


Noticing that the right-hand side is positive and so is the left- 
hand side, we have (12) and (13) after some calculations. 

If e = 5 = 0, we have that - M*l/^ by (16). Noticing 
this relation and applying the same analysis, with (37), we 
have (11)-(13) for this case also. 

(Sufficiency) We employ the control scheme where crk+i 
is determined so that the equality holds in (6) and Uk - 
M*(J/^^j). With this control, we have the equality in (23). 
Moreover, let the quantizer be the optimal one where 
the boundaries are in (15). Then, equality holds in 

(32) also. From (12) and (13), we have that > 0 and (36), 
and then (38). In addition, with (11), the key inequality (37) 
in the necessity part follows and thus, E[cr^] ^ 0. Since 
cr,(./2 > supy^j/^ lyl by definition of the quantizer, we have 
that E[supygj/^ lyp] —> 0 as k ^ oo. □ 


4 General order plants case 


where wf'’ and are given by 


wP 2{\a*\ + edhi. 


wP max |(|a*|-H e,)(l + 

' ie|0,...,riV/21-l)( \ 1^*1/ 


-(Kl-e,)(l-^^|/., 


(42) 


(43) 


We will see in the proof later that iv, is used to express an 
upper bound of the expansion rate of the quantization cells 
enlarged by the parameter Moreover, let the matrix Hr^ 
containing Oi^k, ■ ■ ■, Sn,k be 


0 1 0 


Hr, 


0 0 1 

^n,k ^n-l,k ' ' ' ^l,k 


(44) 


In this section, we consider general order plants in (1), and 
present a control scheme to make the feedback system MSS 
along with a sufficient condition providing a stability test. 
In the course, we will see that results of Markov jump linear 
systems [28] play a central role. A related approach can be 
found in [12], where the known plants case has been studied. 

Given a data rate R = log A and a quantizer {hppp'^, we set 
the control law as follows: In the encoder and the decoder, 
the scaling parameter cr^ is determined by 


Here, is the random vector defined as T^^ := [jk-n+i Jk-n+i 
■ ■ ■ yk\, which is determined by the results of the past n com¬ 
munications. The transition probability matrix P 6 
for the random process (r,t],t is given by, for n>2. 


fn-l 


1 - ^ ^ 0 0 

fn-l 


0 

0 

1 


where (g) is the Kronecker product, and when n - I, 


o-k ^ P(^k) + 


(39) 


P := 


l-q q 
P 1 -F 


at each time k. Here c(-) is the midpoint of an interval. 

Furthermore, the control input is given as Finally, we define the matrix F using //pj and P by 


Myt = - 




(40) 


Next, we introduce some notation. For i - 1,2,...,«, let the 
random variables Oi k be given by 


F :- F\F2, 


(45) 


where F\ P^ ® I„ 2 , F 2 := diag(//p(i) (g) Hpn ,..., //p( 2 ») g) 
Hyu")), and F0\ ..., represent all possible vectors of 
indexed arbitrarily. Let p(-) be the spectral radius of a matrix. 

We are ready to present the main result of the subsection. 


\a*\ + et if yk-i+\ 
Wi if yk-i+i 


0 , 

1 . 


(41) 


Theorem 8. Given the data rate R - log N, the loss proba¬ 
bility p e [0,1), and the quantizer {hppp^, if the matrix F 
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in (45) satisfies 


the control law (39) and (40). This is done by substituting 
(40) into (1) and by the definition (5) of 5/k+i to obtain 


P(F) < 1, (46) 

then under the control law using (39) and (40), the feedback 
system is MSS. 

For the scalar plants case, the inequality (46) holds if and 
only if vwi < 1 and {(1 - q)(\a*\ + ef + (1 - p)w^^ /2 < 1. 
Notice that the first inequality is equal to (37) in the proof 
of Theorem 4 when the quantizer is the optimal one 
Thus, for this case, the condition (46) is tight in the sense 
that if /? e N, /?, g 6 (0,1), and A satisfy (11)-(13), then 
there always exist an encoder and a controller such that (46) 
holds. 

To establish the theorem, consider the Markov jump system 

Zk+\ = Hr^Zk, zo := [o--„+i o-_„ ■ • ■ cro]^. (47) 

The stability of the feedback system can be reduced to that 
of this system. This is shown in the following lemma, whose 
proof is given in the Appendix. 

Lemma 9. If the Markov jump system (47) is stable in the 
sense that converges to the zero matrix, the original 

feedback system is MSS with the control law (39) and (40). 

Remark 10. In [20], a necessary condition for the multi¬ 
dimensional plants case, which is similar to the scalar case 
result, is given. This result has been shown under the setup 
that the structures of state estimators and controllers are con¬ 
strained as (5) and (7). As we discussed in Remark 3, this 
may cause some conservativeness. Therefore, in this paper, 
we do not present the result corresponding to Theorem 4. 

5 Conclusion 




ai,kyk 


' ^n,kyk-n+l bk 




1 

< - 

" 2 




0~ifc+l 

2 


Next, to establish that the stability of (47) implies that [crk}k 
is MSS, we prove the following relation: 

o-k < izk)n for = 0,1,..., (A.l) 


where (•)« is the nth element of a vector. The scaling param¬ 
eter (39) can be decomposed as 




- y n(y[,J/,_,.+i) + |?:W^i)l 

U I 


< y }yP(J{.i3fk-i+\) + 


-i+lj 


(A.2) 


Here, the equality follows from the Brunn-Minkowski theo¬ 
rem and the inequality from applying the triangle inequality 
to the second term. Next, we explicitly compute the width 
p(I7{.i3fk-i+\) of the product of intervals Jli and ffk-M- Re¬ 
call Jli - [a* - e,, a* + e,]. Based on basic results for interval 
products [26], we can obtain 


p (ITliOfk-i+l) - \ 


(|fl*| + e,)pi3/k-M) 

if 3/k-i+i 3 0, 

\a*\p(3/k-i+ \) + ei\3fk-i+\+^k-i+ 1 1 
if 3fk-i+\ t 0 and Ihi t 0, 

2e,max||,y^_,+i|, Wk-i+\\ 

if ^k-i+\ t 0 and HKi 3 0, 


(A.3) 


In this paper, we have addressed a stabilization problem of 
parametrically uncertain plants over data rate limited chan¬ 
nels subject to random data losses. The result for the scalar 
case establishes limitations and trade-off relationships for 
stability among the data rate, the transition probabilities of 
the channel states, and the uncertainty bounds. As mentioned 
in the Introduction, uncertain systems in networked settings 
have not been studied much. We plan to extend our research 
in this area in the future. 

Acknowledgment. The authors would like to thank M. Fujita, 
S. Hara, and R. Tempo for helpful discussions on this work. 
We are also grateful to the anonymous reviewers for their 
comments that helped us improve the paper quality. 


for i - 1,2,..., n. Similarly, the absolute value of the center 
of the interval can be computed as 

\c(J[i}fk-i^l)\ 

0 if J/^,+1 3 0, 

5 (|a*| H- e;)max||.y;i ~i+l I, Wk -i+ll| 

+(|a*l - fi)min||,y^,+i|, 

if ifk-i+i ^ 0 and IHi t 0, 
|a*|max||,y*-/+i|, Wk-i+x'^ 

if }/k-i+\ i 0 and IRi 3 0. 


A Proof of Theorem 8 


We use (A.3) and (A.4) to obtain an upper bound on the ith 
term in (A.2) over all possible 3fk~i+\ as 


Proof of Lemma 9 We first verify that the mean square sta- p{IRi3fk-ivi) + \c(IRi3fk-M)\ < 0i,kO-k-i+i. (A.5) 

bility of [cr^)<: implies that the feedback system is MSS under \b*\ 
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This is shown by examining the following three cases de¬ 
pending on the packet loss process jk-i+i and N. 

(i) jk-i+i - 0: In this case, by construction in Sec¬ 
tion 2 , ^k-i+\ is the entire input range, that is, = 

[-cr^_,+ 1 / 2 ,cr^_,+ 1 / 2 ]. This interval contains the origin as 
an interior point. Thus, by (A.3) and (A.4), we have 

26 

li{Jii}fk-i+\) + k-i+i)\ = (|fl*l + edcTk-i+i. 

\b*\ 


given in (44). Thus, using the nth state of the Markov jump 
system (47), we obtain (A.l). This implies that if 
goes to the zero matrix as k —> 00 , then E[cr^] —» 0 . □ 

From the stability result [28, Theorem 3.9] for Markov jump 
systems, the inequality (46) is equivalent to that (47) is MSS. 
Therefore, from Lemma 9, this inequality (46) is a sufficient 
condition for mean square stability of the original feedback 
system. This concludes the proof of Theorem 8 . 


From (41), we have \a*\ + e, = O^k for this case and hence, 
(A.5) holds. 

(ii) jk-i+i — 1 and N is odd: We must deal with two cases. 
(ii-1) 2^k-i+\ < 0 < 2^k-i+\'- This implies that }fk-i+i is the 

center quantization cell, i.e., its boundaries (J/it-i+i, JZ/t-i+i) 
are {-hi<Jk-i+\,h\(Jk-i+\) (see (4) for the definition of the 
cells). Thus, from (A.3) and (A.4), we have 

n k-i+\) + — |c(J?l,J/i:-,+ l)| = wf^CTk-i+l, (A. 6 ) 

\b*\ 

where is defined in (42). 
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